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Abstract. 

We use a method recently introduced in Barcelo et al, 10.1103/Phys- 
RevD. 83. 041501, to analyse Hawking radiation in a Schwarzschild black hole as per- 
ceived by different observers in the system. The method is based on the introduction 
of an "effective temperature" function that varies with time. First we introduce a non- 
stationary vacuum state for a quantum scalar field, which interpolates between the 
Boulware vacuum state at early times and the Unruh vacuum state at late times. In 
this way we mimic the process of switching on Hawking radiation in realistic collapse 
scenarios. Then, we analyse this vacuum state from the perspective of static observers 
at different radial positions, observers undergoing a free-fall trajectory from infinity, 
and observers standing at rest at a radial distance and then released to fall freely to- 
wards the horizon. The physical image that emerges from these analyses is rather rich 
and compelling. Among many other results, we find that generic freely-falling observes 
do not perceive vacuum when crossing the horizon, but an effective temperature a few 
times larger than the one that they perceived when started to free-fall. We explain 
this phenomenon as due to a diverging Doppler effect at horizon crossing. 
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1. Introduction 

Arguably, the discovery that black holes should evaporate by emitting a Planckian 
spectrum of particles, is the most important result of combining general relativistic 
with quantum mechanical notions [HE]- As it is by now well known. Hawking radiation 
is a kinematic effect (that is, does not directly depend on Einstein's equations), and so, 
much more general than its incarnation in General Relativity. In any phenomenon 
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in nature describable in terms of quantum excitations propagating in a Lorentzian 
geometry with black-hole-like properties, one expects to have Hawking particle (or 
quasi-particle) production. In [3], for example, you can find a catalogue of systems and 
situations outside the General Relativity realm in which Hawking radiation is expected 
to appear. Here, by black-hole-like properties we mean either the presence of trapping 
horizons or just the asymptotic approach towards configurations with horizons [H El |6]. 

When the Lorentzian geometry acquires the form of a stationary black hole, the 
temperature of Hawking radiation becomes constant and proportional to the surface 
gravity of the black hole at the horizon |i2|. Moreover, this temperature only depends 
on this constant, leaving no trace of the history of the black hole geometry formation 
or preparation [2|, [7j. In the General Relativity realm that we shall adopt in this paper, 
this surface gravity and so Hawking temperature, happens to be inversely proportional 
to the mass of the black hole. 

When investigating Hawking radiation, it is customary to work on a simple static 
background geometry, namely the Schwarzschild solution, neglecting any back-reaction 
on the geometry. In order to have Hawking radiation at infinity, it is necessary to 
set the quantum field in a particular quantum state: the Unruh vacuum state [SI E]- 
This state is commonly described as being a vacuum state for observers freely falling 
at the event horizon of the black hole. It is a stationary state, describing a black hole 
which has always been and will ever be emitting radiation at Hawking temperature 
at all times. Therefore, it does not capture the process of the black hole formation in 
General Relativity, in which Hawking radiation switches on in the last stages of collapse. 
Instead, in this paper we are going to analyse the characteristics of a different vacuum 
state, which might be called the collapse vacuum. This vacuum is appropriately chosen 
so as to mimic the process of switching on of Hawking radiation. Initially, there will 
be no radiation at infinity, but at some particular time, some radiation will start to be 
noticeable. After some rather quick transient regime, this radiation will become thermal 
with Hawking temperature (as in this paper we are not going to consider the effects of 
the evaporation, it will not be necessary to distinguish between the terms thermal and 
Planckian, and we will use them indistinctively). Our collapse vacuum interpolates 
between the Boulware vacuum state ^lOj at early times and the Unruh vacuum state at 
late times. 

It is well known that "presence of particles" is an observer-dependent notion [H [11] . 
In this paper we will be specifically interested in analysing how the collapse vacuum state 
is perceived by different observers. In particular, we are going to consider (i) static 
observers at different radial positions of spacetime, (ii) observers undergoing a free-fall 
trajectory from infinity, passing through a fixed radial position at different times, and 
(iii) observers standing at rest at a radial distance until a moment at which they are 
released to fall freely towards the hole. A discussion about the different perceptions 
of the radiation from a black hole in the Unruh state by freely-falling and stationary 
observers can be found in [12]. 

In order to analyse this particle perception, we could use Bogoliubov 
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coefficients [TT], but these, except for a few particular simple cases, are difficult to 
calculate. Instead, we could use the functional Schrodinger formalism as in |T2], but 
it appears difficult to manipulate for non-stationary vacuum states. In this paper we 
will adopt a different strategy. Following [13], we shall define and calculate a time 
dependent effective temperature for each observer along its world line. Then, abusing of 
the language we will describe the perception of the observer as experiencing a thermal 
radiation wind with a time dependent temperature. This effective temperature will 
only be a real temperature in situations in which an additional adiabatic condition 
is satisfied. We will also calculate when this condition is satisfied and when it is not. 
Even in the cases in which the adiabatic condition is violated, and therefore the radiation 
spectrum differs from the Planckian shape, it will be argued that the calculated effective 
temperature is still useful as an estimator of the amount of particles encountered by 
the observer at each instant of time. But it is important to remark that in these 
cases our method cannot give information about the characteristics of the spectrum. In 
any case one should keep in mind that a precise calculation of the particle perception 
would require the construction of a particle detector, taking into account its spatial and 
temporal resolutions. This is, however, beyond the scope of this paper. Nonetheless, we 
will show that the physical image that emerges from our description is rather rich and 
compelling. For instance, it will allow us to clearly see that freely-falling observers at 
the horizon will not perceive vacuum unless (a) they closely follow the collapse or (b) 
they have instantaneously zero radial velocity at the horizon. We will also show that 
this result is not an artifact of the chosen vacuum, but that it is precisely a characteristic 
of the final Unruh state. 

The paper is structured as follows. In section |2| we will start recalling the basic 
background ingredients needed for our analysis. In section |3j we will describe and fix 
the vacuum state for the quantum radiation field of our problem. Then, in section |4| 
after describing what we mean by effective temperature and adiabatic control functions, 
we will study in detail the radiation perception associated with the already enumerated 
different observers. We will devote section [5] to the explanation of one of our most 
important results: the non-zero perception of radiation by freely-falling observers at the 
horizon, against the opposite naive expectation. After a brief comparison of the results 
with those that one would obtain in a pure Unruh state (section |6]), we will finally 
summarize and conclude. 

2. Preliminaries 

2.1. Geometrical set up 

We shall restrict ourselves to the analysis of radiation emission in a Schwarzschild black 
hole. In units where G = c = 1, the Schwarzschild metric reads 




(1) 
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with m being the mass of the black hole. Moreover, we will only use the {t, r) sector 
of the geometry. In double null coordinates, the corresponding (1 + l)-metric can be 
written as 



ds' ^-y-- —j d^id^;. (2) 
Here, 

u:^t-r*, v:^t + r*, (3) 

with r* being the "tortoise" coordinate 

r* :=r + 2mlnf— -iV (4) 
\2m / 

2.2. Quantum field theory 

We are going to analyse the behaviour of a quantum field in this spacctime. For 
simplicity, we will only consider a Klein-Gordon massless scalar field (analyses with 
other fields would yield qualitatively similar results). In this paper, we will only consider 
spherically symmetric configurations and ignore any backscattering of the scalar field 
in the geometry, so that our analysis will be equivalent to that in a (1 + 1) -dimensional 
conformally invariant theory (here on, wc will stick to a pure 1-1-1 analysis). Now, the 
1-1-1 Klein- Gordon equation in null coordinates reads 

so that its general solution can be written as 

cl> = f{u) + g{v). (6) 

Owing to the conformal invariance of the 1 + 1 Klein-Gordon equation, any relabelling 
= Pvs{U), V = qvsiV) (the subscript 'vs' in these functions stands for vacuum state; 
the reason for this notation will be understood later on) can be associated with a mode 
decomposition of the quantum field of the form 

= du' [du^, + di<t>i + ai}{cr + di'}{<t>i,r\ , (7) 

where the normalized modes in this expression are 

This can easily be seen by using the following form of the Klein-Gordon scalar product 
(01, 02) -i j dU<P^ du <t>l + j dl^0i dv ■ (9) 

Once the functions u — Pvs(^), v = q'vs(^) have been defined, one can select a vacuum 
state |0') by requiring that it satisfies 

a^,|0') = 0, a^,|0') = 0. (10) 
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For any other choice u = poh{u), v = qoh{v), we can perform an analogue mode 
decomposition and also define a vacuum state, 

am=0, a:\0)=0. (11) 

In general, the vacuum states |0') and |0) will be different. Thus, an observer for which 
the vacuum is |0) might detect particles corresponding to the vacuum |0'). 



2.3. Bogoliubov coefficients 

The particle content of the vacuum |0') as seen by observers with an unprimed vacuum 
notion can be calculated through Bogoliubov coefficients. In this paper we will only 
worry about outgoing radiation, so we will forget about the 0^ modes, which correspond 
to rays going into the black hole. For each frequency cj, we have (see [TT] ) 



(0'|iV::|0') = {Q'\a:}al\Q') = j Au' , (12) 

where is the Bogoliubov coefficient defined by the scalar product 

Pu.':=-{<PlA€r). (13) 

with 0^ defined analogously to ([s]). Plugging the expressions for these modes into ([9]), 
one arrives at 

(3^^, = ^ (uj I dMe-'('^'^+'^'^) -u' I df/e-^('^'^+^"M , (14) 

A-nyJojuj' \ J J J 

or, integrating by parts the second term, at 

P..' = ^\l^, Uue---'^^%-\ (15) 



27r V W _ 

(where we have eliminated an irrelevant boundary term). Therefore, to determine the 
particle content, one only needs to know the relation U{u). 



3. Vacuum state choice 



We are interested in analysing how Hawking radiation is perceived by different observers 
in spacetime. In order to have Hawking radiation in the first place, we need to select an 
appropriate vacuum state. Instead of making the obvious choice of the Unruh vacuum, 
we shall choose as vacuum state a non-stationary state which somehow interpolates 
between the Boulware state (with no radiation at both past and future asymptotic 
infinities) at early times and the Unruh state at late times. This dynamical (or non- 
stationary) state in a Schwarzschild static background will mimic what happens in 
a collapse process, in which the black hole is generated from an initially (almost) 
Minkowskian spacetime. As opposed to what happens in the Unruh vacuum, in this 
situation Hawking radiation will not be present at early times, but will be switched on 
at a particular ignition time (this would correspond to the time of the formation of the 
black hole). 
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To impose this vacuum state to the quantum field, we will consider a time-like 
geodesic observer falling into the Schwarzschild black hole from infinity. We will 
then calculate how this observer would label in a natural way the different outgoing 
{u = const) rays that he encounters in its way towards the horizon. This new label 



U 



P. 



-}{u 



will then be used to define the vacuum state via a positive-frequency mode 



expansion, as explained in section 2.2 [see ([7j)]. 

So let us first calculate the radial time-like geodesic trajectories by solving the 
geodesic equation 

/ j„ \ 2 / \ 2 



rr 



rr 



d r 

With this aim we will take into account the relation 

2 / o \ -1 



dr \ ■ 

d7y 



dt Y 



2m Y 
r ) 



1+1 



2m Y 
r / 



dr Y 
d7 J 



(16) 



(17) 



that follows from the spacetime interval (equivalently, we could use directly a variational 
principle for the spacetime interval). Here, r is the proper time of the trajectory, and 
r^j, and are the only non-vanishing Christoffel symbols relevant for the calculation 
of radial trajectories. 



m / ^ 'im \ ^ m 



^ 2m\ 
r J 



r J 



Replacing them in the radial geodesic equation yields 
d^r m 

d7"~^' 

which has a form resembling Newton's gravitational force. Integrating once, we obtain 

1/2 

1 - - 1 , (20) 



(19) 



2m 

r 



r 



1 



2m Y 
r / 



1 



2m \ 



1/2 



(21) 



dr 
d7 

dt 

dr \^ T ) \ tq J 

in which we use the radius tq at which the velocity vanishes as an integration constant. 
For the time being, we are interested in describing a free-fall trajectory starting at radial 
infinity, that is, with tq oo. Then, integrating the radial equation we obtain 

2/3 



3V2m 



(ro - r) 



(22) 



with To being an integration constant. Note that r runs from — oo to Tq, at which time 
the trajectory reaches r = 0. Before that, at tq — 4m/3 the event horizon is crossed. We 
can also find the behaviour of the time coordinate t with respect to r. Indeed, dividing 



equations (20) and (21) (with tq oo), we obtain 

dr /2m 2m' 

dt ~ ~V~ 



1 - 



(23) 
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which can be integrated to yield 
t = to — Am 



r 



1 / r \3/2 1 

+ - + - In 

2m 3 \2mJ 2 



^^/rJ(2m) 



v/r/(2m) + 1 



(24) 



where to is an integration constant. Finally, replacing the relation (22) for r(r) in this 
expression, we obtain 



t = to — 4m 



3(ro - r) 
4m 



1/3 



+ ^ + -ln 
4m 2 



(3(ro - r)/4m)i/3 - 1 
(3(ro - r)/4m)i/3 + 1 



(25) 



Now, having the pair (t, r) as a function of r in (22) and (25), we can also write the 
quantity m in (|3| as a function of r, 

\ 1 2/3 



U 



Am <| 


'3(ro-r)" 


1/3 ^ 


'3(t„-t)' 




4m 


+ 2 


4m 



+ 



To 



4m 



+ ln 



3(ro 



4m 



1/3 



- 1 



(26) 



An observer following this world line can naturally use its proper time r to label 
the different rays he encounters in its way towards the horizon, namely, he can make 
the assignment of the label U (u) = r, given by the inverse of (26), to the ray u = const 
that hits the observer at proper time r. Then, defining Un := tq — 4m/3 and choosing 
the two constants to = so that the two labels U and u are synchronized at u — — oo, 
in the sense that U = U{u — )■ — oo) — )• -u, the final relation between both labels is 



u = p^^{U) := Un + 



+ 




4m 



(f/H - f/) + 1 



(t/H -U) + l 



2/3 



1/3 



1 



(27) 



Note that, while the range of the original null label is m G (— oo,+oo), the range of 
the new null label is f/ G {—oo, Uu), where Uu marks the moment at which the falling 
observer crosses the event horizon. As explained before, associated with this choice 
U 

subject of our analysis 



p^^{u), we have a vacuum state, and it is this vacuum state which will be the 



4. Radiation perception for different observers 

In this section, we will investigate what is the particle content of the vacuum state 
defined above for different types of observers, for which the notion of vacuum might or 
might not coincide with this one. More specifically, we are going to define yet another 
label u = p~b (m) associated, in each case, with a specific type of observer. Once we have 
characterized the observer by the null-ray labelling function, we will find the function 
U = Pvs[Poh{u)] =: U{u) which will inform us about the (observer-dependent) particle 
content of the chosen vacuum state. 
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Variable temperature-like estimator 



As an alternative to the calculation of the exact Bogoliubov coefficients using (15), we 
can estimate the amount of particle content by calculating the function 

:= -— ( — 1 . (28) 



As described thoroughly in [U [T3j, when this function is constant k{u) ~ =: k^, 

over a sufficiently large interval around a given u^^ one can assure that during this same 
interval the system is producing a Hawking flux of particles with a temperature 

fcBT=^. (29) 

If the variation of k[u) is slow, then one can describe the system as a thermal 
emitter with a slowly varying temperature. Whether k{u) varies slowly or not in the 
surroundings of u*, is controlled by an adiabatic condition which, under mild technical 
assumptions [13], reads 



1 


d/t 






dw 





< 1. 



(30) 



However, even in the case in which the adiabatic condition is not satisfied, the 
existence of a non-zero k[u) is an indication that there is particle emission (in this case 
the Bogoliubov /3 coefficients will not be zero). The spectrum of this particle content 
does not follow a precise Planckian proffie anymore, but still we will take k{u) as an 
estimator of the amount of particles perceived by the observer. In this work, for short 
we will refer to «:(«) [or more precisely to n,{u) / (27r)] as a variable temperature, although 
strictly speaking this term is correct only in the intervals in which the adiabatic condition 
is satisfied. We will call 



1 



du 



(31) 



the adiabatic control function 



4-2. Static observer at a fixed radius 

Let us start by analysing the particle content as seen by a static observer sitting at 
a radius r^. In order to find a new labelling u = Pso{u), now associated with this 
static observer, we will repeat the same steps that we did before for the definition of the 
vacuum state. (Here on we will replace the generic subscript 'ob' by a subscript denoting 
the particular type of observer; in this case 'so' stands for static observer.) The world 
line of an observer keeping its radial position can be described as r(r) = = const and 

t=(l-^-^y'^\r-r,), (32) 

so that 

/ 2m\"^/^ 

u=\l-—j (r-ro)-r*(r3). (33) 
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As before, an observer following this world line can naturally use its proper time r to 
label the different rays he encounters while standing in its radial position, namely, he can 



make the assignment of the label u = p\ 



U 



T given by the inverse of this equation 



to the ray u = const that hits the observer at proper time r. Upon synchronization, we 
obtain the relation between both labels 

-1/2 



U = Psoiu) 



1 



2m 



u. 



(34) 



Note that here one cannot demand m — ?• m in the past infinity, as we did in section [3j 
The synchronization in this case reduces to eliminate an irrelevant additive constant. 

At this point, we can already construct the relation U = U{u) we were after from 
the equation u = Pvs{U) = Pso{u). The complexity of (27) does not allow us to write 
this relation explicitly, but we can write its inverse. Using (27) and (34), we obtain 

4m 



u 




— Am 



4m 

2/3 



(f/H - f/) + 1 + 



(Un - f/) + 1 



1/3 



(f/H - f/) + 1 



1/3 




(f/H - f/) + 1 



(35) 



One can study the behaviour of this function at the past and future infinities, m — )■ — oo 
and u — )■ +00, respectively. In the asymptotic past, the linear term on the right-hand 
side is the most important, and this relation becomes 



U 



1 



2m Y 



-1/2 



U, 



when u — )■ — oo 



(f/^ 



-oo 



(36) 



By looking at the definition of k{u) in (28), one can easily see that in this regime k = 0, 
so that there is not radiation at all. 

On the other hand, at future infinity the behaviour of U{u) is dominated by the 
logarithmic term in (35), yielding an approximate expression 

2m\-^/^ 



f/^f/H 



4m 

IT 



exp ( - — (f/n — 6m) ) exp 



\4m 



1 - 



1 



4m 



(f/H - 6m) 



exp 



rs J 



1 



+ 1 



2m \ 



-1/2 



U 

4m 



~ f/n — 4m exp 

when u ^ oo iJJ ^ f/n). 
From this expression, it is easy to recognize the well known asymptotic form |2] 

U = Un- v4exp(-Kau), 

(where Uii,A,Kg, are constants), which is associated with a Planckian spectrum of 
particles with temperature k^T = fi;a/(2vr). Thus, when the static observer waits long 



(37) 
(38) 
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enough, it perceives a Hawking flux with temperature 



Ts J 



-1/2 



(39) 



This is precisely the asymptotic Hawking temperature k^T^i = KH/(2vr) = l/(87rm), 
multiplied by a gravitational blue-shift factor which runs from unity for observers close 
to infinity, to arbitrarily large values for observers standing at positions closer and closer 
to the event horizon. 



4-2.1. Numerical k,{u) Let us describe now the particle perception of one of these 
observers throughout its entire life outside the horizon. From the definition of k{u) 



in (28), it is easy to see that we can also write 



d^u / du Y 



U(u) 



Then, evaluating the derivatives we can obtain a close form for n{u) 



K[U) 



1 

4m 



-4/3 



(40) 



(41) 



Solving numerically the relation U{u) from (35), we can plot k{u) for different static 
observers (figure [T]). 




Figure 1. Temperature as a function of u for different static observers at r-g — {3m, 

4m, 10m, 20m) [curves depicted respectively as ( , , — • — , )]. We use 

2m = 1 units. 

In these graphs, it is clearly seen how the vacuum state is such that in the past 
it does not contain particles, but at some instant of time it starts to heat up reaching 
asymptotically a final temperature. As we mentioned at the beginning, the vacuum 
state that we have chosen in a static background mimics what would happen in the 
formation of a black hole through gravitational collapse. Only once the black hole is 



very close to formation one enters in the asymptotic regime (37), and Hawking radiation 



is switched on. We can estimate when this happen, but for that let us first describe the 



behaviour of the adiabatic control function in (31). 
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4-2.2. Validity of the adiabatic approximation From the definition of e in (31 ), we have 

1/3 



e(u 



dU 



, du 
dU 



-1 



4 



U{u) 



4m 



(Uu - Uiu)) + 1 



1 



(42) 



Thus, e decays from infinity at the asymptotic past, to zero when the observer that fixes 
the vacuum reaches the event horizon. This late time behaviour is consistent with the 
fact that then the quantum field is finally switched on, so that the radiation that a static 
observer sees is perfectly thermal. On the other hand, in the asymptotic past, e becomes 
unbounded, but this does not mean a failure of the adiabatic condition, because k — ?■ 
there. Again, we can numerically plot e{u) for different static positions (figure |2|. We 
can see that its value crushes to nearly zero when the temperature enters its asymptotic 
regime. 



10 - 




Figure 2. Adiabatic control function e as a function of u for different static observers 

at Ts = (3to, 4m, IOto, 20to) [curves depicted respectively as ( , , — • — , 

)]. We use 2m = 1 units. 



Now we can introduce two different referential times. On the one hand, we can 
define an ignition time of the black hole as the moment at which k = 0. If one 
calculates this ignition time in terms of U, it does not depend on the value of r^. This 
is perfectly reasonable, as the label U corresponds to the single free-fall trajectory used 
to select the vacuum state. It is not difficult to see that this ignition time corresponds 
to t/jg = Uh — (676/1029)m, or equivalently to the time when the free-fall trajectory 
reaches Hg ^ 2m + 0.61m. On the other hand, we can define a thermality time as the 
moment from which the outgoing radiation can be assumed to have a Planckian shape. 
This thermality time can be calculated in terms of U through the condition e{Uth) = ef, 
meaning that the adiabatic control function has reached a fiduciary value ef, that here 
we will take equal to 0.01. This specific condition happens at Uth — Uu — 0.01m, 
which corresponds to the time when the free-fall trajectory reaches rth — 2m -|- 0.01m. 
Let us give you some figures. If a neutron star with 1.5 times the mass of the Sun 
(Schwarzschild radius equal to 4.43 km) and a radius of 12 km started to undergo a 
free-fall collapse at Schwarzschild time tini, it would reach rjg ~ 2m -|- 1.357 km at a time 
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tig ~ tini + 0.13 X 10"^ s, and rth ~ 2m + 0.022 km at a time tth - ^ini + 0.20 x 10"^ s. 
In this way, one can directly see that the detection of Hawking quanta at infinity does 
not involve any long delay associated with the freezing of the collapsing structure as it 
approaches the horizon formation. 

Finally, let us comment that for different observers the previous two times are seen 
to happen at different instants of their proper times (their respective u labels). The 
closer to the horizon, the later the observer will see Hawking radiation to switch on and 
to become Planckian. 



Jf-.S. Observers freely falling from infinity 



In this subsection, we will analyse particle creation as seen by observers which are 
freely falling from infinity. To find the relevant re-labelling u = pmo{u) (where 'ffio' 
stands for falling-from-infinity observer), we have to use the form of the radial time- 
like geodesic that starts at spatial infinity with zero velocity. However, recall that the 
calculation done to fix the vacuum state of the field in section [3] was precisely using these 
time-like trajectories. In that calculation, we had two constants of integration to play 
with, namely to and tq. We will now consider trajectories differing from the reference 
trajectory in section |3] by a temporal delay Ato. Instead of using t^ = tq = Un + 4m/3, 
we will use new values to = = + 4m/3 + Ato = + 4m/3, where un := Un + Ato. 
We then define the new relation u 



u = pmo[u 
which explicitly reads 
4m 



u 



un + 



— 4m 



Psioiu) by [see ( [27| )] 
{u - Ato) + Ato, 



(43) 



4m 



[un -u) + l 



1/3 



+ 



-1 2/3 



4m 



Am 



(uh - m) + 1 



+ In 



4m 



u] 



'uh -u) + l 



1/3 



(44) 



Now, as we did for the static observer in section |4.2[ we compare both labellings, U 

(45) 



in (27) and u in (44), and obtain an implicit relation U{u) from the equation 

P^s{U) =Pmo{u), 
which again explicitly reads 



— 4m 



4m 



(f/H -U) + l 



1/3 



+ 



4m 



{Uu - U) 



1 2/3 



4m 



(f/H -U) + l 



In 



Ato - 4m 



4m 



(mh -u) + 1 



4m 

nl/3 



(f/H -U) + l 



1/3 



4m 



(mh -u) + 1 



2/3 



+ - 



4m 



(uh -u) + 1 



+ In 



4m 



(uh -u) + 1 



1/3 



(46) 
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Solutions to this non-algebraic equation must be found numerically. Nevertheless, 
we can examine the past infinity and event horizon limits analytically. In the past 
infinity, we obtain U ^ u, owing to the synchronization we imposed. Near the event 
horizon, the logarithmic functions provide the leading contribution and therefore, 

3 



4m 



-Ato/{4m) 



4m 



[uu - u 



1 



1/3 



1 



1 



(47) 



The function k{u) can be readily obtained from (43) and (45). Indeed, taking the 



derivative of (|45|) with respect to u, we obtain 



dU 
du 



Ato) 



du 



dU 



Uiu) 



and differentiating once more 



d^U 
du^ 



du^ dU^ 
so that k{u) can be expressed as 



(-Y 

u{u) \ duj 



dpUU) 



dU 



U{u) 



dU' 



U{u) 



duy dv 

du J 



[u 



Atn 



du' 



dp^sju - Ato] 
du 



-1 



(4J 



(49) 



(50) 



Therefore, the final result is 



k{u) 



4m 



n -1 



4m 



X 



(uh -u) + 1 



(3/4m)(MH 



4m 
u) + l 



(mh -u) + 1 



nl/3 



- 1 



-1 4/3 



_(3/4m)(f/H-f/(M)) + l. ^^^^ 

Note that this exact expression is explicit if we considered it as a function of both 
u and U (although this does not mean that it is a function of two variables). The 
dependence on Ato is hidden in the implicit relation between both variables. We can 



use the limiting expression of U{u) at the event horizon (47) (which is exact there) to 



find what temperature would the freely-falling observer characterized by the delay Ato 
see when he crosses the horizon. This temperature happens to be 

1 



Knc{AtQ 



m 



(1 



-Ato /(4m) 



(52) 



Therefore, the perceived radiation will run from an initial zero temperature to nearly 
l/(27rm), which corresponds to four times the standard Hawking temperature, provided 
that the waiting time Ato is several times bigger than 4m. This is an interesting and, 
at first, puzzling result. How can one interpret it? In section |5| we will offer and 
explanation of this result, but in order to do that we still need to introduce several 
additional material. 
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4-3.1. Numerical k,{u) We can numerically plot the entire relation k{u). Results 
are shown in figure [s] For small values of Atg (less than approximately 20m) the 
temperature perceived by the freely-falling observer is nearly zero during most of its 
trajectory towards the horizon. Only when approaching the horizon, the temperature 



increases to reach a peak value given by (52). For large enough values of Atg, before this 



final increase in the temperature, there appears an intermediate plateau in which the 
temperature is nearly constant. It is easy to understand this behaviour: if the freely- 
falling observer is still far from the horizon when Hawking radiation is switched on, at 
some moment it will start detecting this radiation, which will stay nearly constant till 
the observer becomes close to the horizon, where additional physics will come into play 
(see section [5]). The effective temperature value along the plateau is somewhat higher 
than the Hawking radiation temperature, and slightly increasing. As we will see, this 
is mostly due to a Doppler blue-shift associated with the radial velocity of the falling 
observer (see also section Isl). 




Figure 3. Temperature as a function of u for different freely-falling observers from 
infinity with delays Ato = (2m, 20m, 100m, 200m) [curves depicted respectively as 
( , , — • — , )]. We use 2m = 1 units and uh = 0. 



It is interesting to look at how good is the functional approximation (47) when 



plugged into (51), as compared with the exact result. This comparison is shown in 



figure |4] for a long delay Ato. By construction, this approximation works perfectly well 
in the surroundings of the horizon crossing. However, it is remarkable that it also fits 
very well the plateau region. 

On the other hand, this approximation is not useful when Ato is not large enough 



see figure 5|. The asymptotic approximation (47) correctly reproduces the horizon- 



crossing value as should be by construction, but at early times it deviates from the 
exact curve rather quickly. 



4.3.2. Adiabatic approximation validity Having found an expression for k in terms of 
u and U (u), we can use it to obtain an expression for the adiabatic control function e{u) 
as defined in (31). The resulting expression, also explicit when considered as a function 
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Figure 4. Temperature as a function of u for a freely-falling observer from infinity with 

A^o = 200r7i using the exact expression (curve depicted ) and using C/near-hor(M) 

approximation (curve depicted ). We use 2m = 1 units and wjj — 0. 
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Figure 5. Temperature as a function of u for a freely-falling observer from infinity 

with Ato = 4™ using the exact expression (curve depicted ) and using C/ncar-hor (w) 

approximation (curve depicted ). We use 2m = 1 units and mh = 0. 



of both u and f/, is rather involved, and we will not include it here. But, as we did 



with K, we can also replace the relation t/near-hor(w) (47) valid in the near- horizon limit 
to find an explicit expression for e{u) valid there; finally, we can take the m — )■ mh limit 
and find its value at horizon-crossing. The result is 



eHc(Ato 



+ 



7 



(53) 



4 [gAto/(4m) _ X] ■ 

This expression runs from 3/8 when Ato is several times larger than 4m, growing to 
infinity when Ato goes to zero. This last value is comprehensible from the definition 
of e, as it is the temperature itself that goes to zero when the freely-falling observer 
approaches the referential one. The adiabatic condition e <^ 1 is never satisfied at 
horizon crossing. However, it is remarkable that, for Ato sufficiently large, it is always 
still smaller than unity, so that the perceived particle spectrum should not be very 
different from a Planckian spectrum. 
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Nonetheless, this is just the horizon crossing value. We can see using the numerically 
evaluated exact result (figure [6]) that, for long waiting times, there is part of the 
trajectory for which the adiabatic condition is with no doubt valid, as there e crushes 
to nearly zero. This of course coincides with the plateau region (in the cases in which it 
appears), in which k takes a nearly constant value. Once the final increase in k shows 
up, the adiabatic condition starts to be violated, although in a rather mild manner. 




-80 -60 -40 -20 

u 

Figure 6. Adiabatic control function e as a function of u for different freely-falling 
observers from infinity with delays Ato = (lOOm, 200r7i) [curves depicted respectively 
as ( — • — , )]. We use 2m — 1 units and uh = 0. 



4- 4- Observers freely falling from a finite radius 



The last situation we shall consider in this paper is the particle perception of a freely- 
falling observer left to fall from a radius ro, with zero initial velocity at some waiting 
time Ato (in Schwarzschild time) after the reference trajectory used to define the vacuum 
state had crossed this radial position. Again, we start by calculating this new trajectory. 



For that, we have to integrate (|20|), but now keeping tq a finite constant. This yields 

1/2 r 



tq-t = -ro 



2m 



1 



+ arctan 



r 



1/2 



(54) 



r / To 

where Tq represents the time at which the observer starts to fall (note that, in all 

expressions in this subsection, 2m < r < ro). The time of horizon-crossing for this 
trajectory is 

' 1/2 rn r / rn \ 1/2' 



ur '■= To + 2rri 



^^0 

2m 



( — 1^ H — arctan 

\2m / 2m 



By means of (20) and (|21|), we obtain 



1/2 



rro_ 
\2m 



'1/2 



(55) 
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which can be readily integrated to yield 




1/2 r 



2m 

ro 



To 

1/2 



+ (1 + 

r 



4m \ 



^0 / 

1/2 



arctan 



To 
r 



1/2 



r 
2m 



2r 



(57) 



where to is again an integration constant, in this case representing the time at which 
the observer starts to fall from tq. If we call t^s{r) the expression in (24), we can define 
to as a function of ro, to aiid the delay time Atg 



to = tvs(ro) + Ato 

= to — 4:m 



ro 



2m 3 \2mJ 




(5J 



Substituting in (57) the function r(r) which is implicitly defined in (54), we can find 
the function t(r). 

Let us now construct the spacetime trajectory of an observer that remains at a 
fixed position r = vo until r = fo, at which time he starts to fall freely towards the black 
hole. This trajectory (t(r),r(r)) is made out of two pieces jointed together at r = tq: 
the first piece, for r < tq, is such that the radius is fixed r = ro and the Schwarzschild 
time satisfies 

2m\-'/' 



tr) 



1 



(^-T^o); 



(59) 



the second piece, for r > fo, is constituted by the previously found implicit expressions 
r{T]ro,fo) in^(|54)) and t(r(r; ro, fo); tq, to, Ato) = t(r; ro, fo, to, Ato) as defined in (|57|; 
the to and Ato dependence of this last expression comes from the dependence of to on 
to and Ato in (58 ). 

As explained in previous sections, substituting r by u, using an appropriate 
synchronization between to and fo, and substituting fo by un as defined in (55), we 
can construct the function 



u 



Pfivo{r) := t(r;ro,MH, Ato) -r*(r), (60) 

where 'ffro' stands for freely- falling-from-a-radius observer. Finally, the U{u) relation 
we were looking for can be found from the implicit relation 

u = Pvs{U) = Psro{riu)). (61) 

Our next step is to find the form of k,{u) for these observers. Taking derivatives of 
(61) with respect to u, we obtain 

d?7 _ / dpvs dpffro dr 
du \ dU ) dr du ' 



(62) 
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where dr/du is nothing but (20) with u instead of r, and differentiating again, 

-1 



d^ro / dr \ ^ dpffro dV _ d^Pvs ^ df/^| ' 



dr vdM/ dr dn dU \du J 



dpvs 
df/ 



(63) 



In this way, /t(-u) can be expressed as 

'dpvsA dVs dU ( dpsTo\ d^pffro dr 



K[U] 



dr\ ^ d^r , ^, 
1 — (64) 



dU J dU"^ du \ dr J dr^ du \du J du^ 
Substituting exphcitly the corresponding expressions and simphfying, it becomes 



K[U] 



k{u) 



for u < uo, 




(65) 



for u > Mo, 



(66) 



where r and U must be understood as functions of u (again, these can only be obtained 
numerically). This function has a finite jump at m = Mq (= ^0)5 the moment at which 
the observer at tq is released from its previously fixed position. 

At early times, k follows the same pattern described in section 4.2, in (41). It 



increases from zero and tries to reach the asymptotic value associated with the observer 
at rest at position tq. However, at some moment Uq the observer is released to fall freely. 
In this moment, k, undergoes a finite jump. The value of k, just after the jump, which 
we will denote ^released, is found by evaluating (66) at r = tq, which gives 

'treleased(^ 0) '^to) = ( 1 



X 



1 

Am 
3 



4m 



2m Y 



(f/H-f/o(ro,Ato)) + l 



-4/3 



-(^) 



(67) 



Here Uq := U{uq) is the value of the label U associated with mq, and depends on the 
parameter Ato- A particular case is Ato = 0, which means that the observer starts to 
fall just when the observer that fixes the vacuum state passes by its side. This imphes 
of course Uq = U{ro), as determined by the trajectory in section [3] [it can be obtained 
by reversing ( [22| , noting that U acts as r now]. Doing this one can check that ^released is 
zero, as one should expect. This is because in this case, at the starting point, between 
the observer that determines the vacuum state and our observer there is no difference 
but in their velocity. If the former one sees nothing, the latter shall see "nothing times 
a red-shift factor" , that is, nothing at all. 
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On the other hand, if the observer waits long enough so that the black hole is nearly 
completely switched on (Unruh state), that is, Atg — oo, which amounts to substituting 



2m\ 



2m\ 



1/2 



(68) 



U by f/n in equation (67), then ^released takes the value 

'^Unmh(''^o) := 'treleased(^o)lAin^oo ^ ~a ( + / I ' / 

The limit tq — t- 2m gives zero, and the limit tq — t- oo gives the usual Hawking 
temperature n = l/(4m), as happens in the Unruh state. This agrees with the fact 
that our state reproduces the Unruh state for late times. 

Another interesting quantity we can calculate is the value of the jump in k, A/treiease, 



at the instant of releasing. We can directly subtract from (65) the expression in (67) 
and find the jump A/treiease, 

1 f ^ 2m\ ^^^/2m\^ / 2m\^^^^m 



Ak 



release 



(ro) 



4m 



1 - 



1 



,2 ■ 



(69) 



It is remarkable that it does not depend on U or, in other words, does not depend on 
the waiting time Atg- This is reasonable, as by switching off the rockets the observer 
is just loosing the radiation that comes from its own acceleration, and this acceleration 
only depends on the radius it was staying at. In fact, in the last expression we can 



directly spot the surface gravity associated with position ro [see (19)] multiplied by 
its corresponding gravitational blue-shift factor. We can conclude then that A^reiease is 
nothing but Unruh radiation associated with an accelerating observer [SlIII]- In figure [7| 
we plot the value of k around the release point for different waiting times. There, one 
can directly see that the jump in the temperature is always the same. 
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Figure 7. Temperature around the releasing point as a function of u for tq 

delays Atg = (0, 12m, 24m, oo) [curves depicted respectively as ( , 

)]. We use 2m — 1 units and ug = 0. 



8m and 



Now, let us study the radiation perception at horizon crossing. The procedure is 



analogue to the one followed in section 4.3 in (61) we identify the diverging terms at 



the event horizon and, from the resulting equation, isolating U we find an expression 
f/near-hor('^) Valid near the event horizon. The expression is quite complex and not really 
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useful by itself, but again, plugging it into (66) we obtain an expression /t(r) valid near 



the event horizon. Finally, performing the limit r — > 2m we find the horizon crossing 
value khc, which in this case is 

1/2 / ^^77^ . \ 1/2 

khc('^o, Ato) 



3/2 




X exp 



3 \2m) 



(4m ~'~ 



1/2 



(70) 



When the observer waits for a sufficiently long time before being released from its 
position (mathematically, Ato — ?• 00), it simphfies to 

1 



Ato- 



(71) 



which reproduces the result k = 1/m for an observer freely falling from infinity [see the 



Atg — 00 limit of (52)]. It is also interesting to analyse the limit of fi;Hc('"0) Ato) when 



ro — )■ 00. The result is 



ro—^co 



m 



(72) 



which, surprisingly at first sight, does not depend on Ato. That is, we cannot reproduce 



the entire formula (52) by taking the rg — ?■ 00 limit. The reason is that the trajectory 



used to fix the vacuum state does never coincide with the trajectory of a released observer 
even if the release point tq is taking close to infinity. Although the velocity difference 
of these two trajectories at tq tends to zero when tq — )■ 00, they arrive at positions 
close to the horizon at very different Schwarzschild times. The freely-falling-from-a- 
radius observer for far enough releasing radius arrives to the horizon always so late with 
respect to the reference trajectory, that it sees the black hole as virtually switched on 
(in its final Unruh state). 

Plotting expression ( 70 ) (figure 8| , we can see that the parameter Ato plays almost 
no role in the temperature seen when crossing the horizon. 



4.4-1- Numerical k{u) Let us now plot some examples of exact numerical resolution 
of k{u) for different trajectories. The aim is just to have a visual and qualitative 
understanding of the global evolution of the temperature, as well as to see how the 
calculated exact limits and behaviours are reproduced. First, we will plot the value of 
K,{u) for observers starting to fall from the radius tq = 10m with different waiting times 
(figure [9]). Graphs show how k{u) starts from very different values when the observer is 
released, but converges to virtually the same final value no matter the waiting time, as 
the starting radius is far enough. 

If the radius tq is quite near the event horizon {tq = 4m in figure 10), there is a 



small but appreciable difference in the value of k{u) when crossing it, depending on the 
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Figure 8. Temperature when crossing the horizon as a function of rg for delays Aip = 

(0, 2to, bra, oo) [curves depicted respectively as ( , , — • — , )]. We use 

2m = 1 units. 
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Figure 9. Temperature as a function of u for tq = lOm and delays A^q = (0, 24m, 

36m, oo) [curves depicted respectively as ( , , — • — , )]. We use 2m = 1 

units and Uq = 0. 



waiting time. We can also see that the jump in when releasing the observer is more 
important in this acceleration to stay at the fixed radius is much bigger here. 

On the other hand, if the radius is really far from the event horizon (ro = 30m 
in figure 11), the jump in n{u) is completely negligible. Different waiting times in this 
case just determine when you start perceiving the radiation, which can happen before 
or after the releasing point. The plateau found in section 4.3, and of course the final 
peak, are again reproduced, as should happen for a large radius. But note that, as we 
already argued, it is not possible to reproduce the cases in section |4.3| for which Ato is 
small. The limit tq — )■ oo here yields the same limit than Atg — ?■ cxd in section 4.3, no 
matter the value of Ato here. 



4--4-2. Validity of the adiabatic approximation Repeating what we have done for the 
observers studied so far, here we can find a formula for the adiabatic condition €{u), in 
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Figure 10. Temperature as a function of u for tq = 4m and delays Ato = (0, 4m, 

10m, oo) [curves depicted respectively as ( , , — • — , )]. We use 2m = 1 

units and Uq = 0. 




Figure 11. Temperature as a function of u for ro = 30m and delays Ato = (0, 

60m, 140m, oo) [cm-vcs depicted respectively as ( , , — • — , )]. We use 

2m = 1 units and Uq = 0. 



this case as an explicit function of both U = U{u) and r = r{ti). Also, once more we 
can plug the expression f/ncar-hor('") valid near the event horizon, and find a function 
e(r) valid there. Finally, we may take the limit r — )■ 2m and obtain the value chc at 
the horizon crossing. Unfortunately, in this case the approximation ?7near-hor(*"), and 
the exact expression e{U,r) are extremely complex in their functional forms. Thus, it 
seems neither possible, nor interesting, to find an explicit expression of enc in terms of 
the parameters ro and A^o- Nonetheless, it is still possible to obtain an expression for 
e(r, ro) (not only in the horizon crossing) when A^o — >■ oo. 
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In this expression, one can now take the hmit r — )■ 2m and obtain enc when Ato 
in terms of Tq, 



Ain 



oo 



(74) 



This resuh runs from 3/8 when tq ^ 2m, to infinity when tq — )■ 2m. With the first 



hmit 3/8, we reproduce again the resuh obtained in section 4.3 for long delays [see the 
limit of (|53|)]. The second divergent limit reflects the fact that observers released from 



positions near the horizon by no means perceive thermal radiation. In figure 12, we plot 



the function in ( 74 ) , together with some numerically solved graphs of enc for zero and 



finite waiting times Ato- Note that the limit 3/8 when rg — )■ oo is common for all the 
graphs. 
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Figure 12. Adiabatic control function e wlien crossing the horizon as a function of 

To for delays A^o = (0, 2to, 5m, oo) [curves depicted respectively as ( , , 

— • — , )]. We use 2m = 1 units. 



Also, we can find an expression for ereieased at the beginning of the falling part of 
the path, again for Ato — ^ C)0. We obtain 

-1 



^released 



Ain 



\2mJ 



(75) 



This yields zero for ro ^ 2m, and infinity when tq — )■ 2m. The zero value appears 
because, at large radius and for large waiting times, what the observer initially sees 
is a perfect Hawking radiation. The divergent value for ro near the horizon has the 
same explanation as in the case of enc- for trajectories starting too close to the horizon, 
the adiabatic approximation is never valid. Neither at the beginning of the fall, nor at 
horizon crossing. 

Finally, we will plot the value of e{u) for the same cases we plotted k{u) in 
section 4.4.1 (figures 13, 14, and 15). As a consequence of the sudden change in k 



at the release time, these graphs should exhibit a Dirac-delta peak in there. This has 
been avoided in the figures for simplicity. Apart from that, graphs need few comments. 



They of course reproduce the limit values already found. As happened in section 4.3 
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when plateaus appear e becomes nearly zero, and so the radiation is perfectly thermal, 
with the temperature determined by k. For earlier parts of the trajectory, e diverges, 
but as we have already explained, this is only reflecting that the temperature there goes 
to zero. 
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Figure 13. Adiabatic control function e as a function of u for vq = 10m and delays 

Ato = (0, 24to, 36to, oo) [curves depicted respectively as ( , , — • — , )]. 

We use 2m = 1 units and uq = 0. 




Figure 14. Adiabatic control function e as a function of u for ro = 4m and delays 

Ato = (0, 4m, 10m, oo) [curves depicted respectively as ( , , — • — , )]. 

We use 2m = 1 units and uq = 0. 
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Figure 15. Adiabatic control function e as a function of u for tq = 30m and delays 

A^o — (0, 60to, 140m, oo) [curves depicted respectively as ( , , — • — , 

)]. We use 2m = 1 units and uq — 0. 



5. Final increase of the effective temperature 

At this stage, we have aU the necessary ingredients to explain the, at first sight puzzling, 
final increase of k for freely-falling observers. For sufficiently long waiting times, we are 
sure the vacuum state is indistinguishable from Unruh state (see also next section). As 
standard lore in the literature, one can read that this state is a vacuum for observers 
freely falling at the horizon. However, we have here several of these observers and 
the value of khc is non-zero for almost all of them: recall for instance the value of 



khc for observers freely falling from infinity ( 52 ) and for observers freely falling from a 



finite radius (70). The only exceptions are the observer closely following the reference 
trajectory and the one with zero instantaneous radial velocity at the horizon. 

Indeed, once Hawking radiation is switched on, there is only one effective 



temperature that does vanish at the horizon: the Unruh temperature Kunruh/(2vr) (68). 
As we explained, this is the k associated with an observer in free fall at a radial position 
r and with an instantaneous zero radial velocity in there (we will call this observer an 
Unruh observer), 

.„,„.(r) = -(l + -j(l--) . (76) 

Let us now make the reasonable hypothesis that all the physics concerning the 
behaviour of k[u) can be understood in terms of local properties of the observer's 
trajectory: its position, its velocity and its acceleration. We will see that this is indeed 
the case. Under this hypothesis, the only characteristic that distinguishes the Unruh 
observer from any of the other freely-falling observers is that the latter have non-zero 
radial velocity. The difference in radial velocities between the freely-falling observers 
and the Unruh observer implies that their respective radiation perceptions have to be 
related by a Doppler shift factor. In fact, the Unruh observer is such that in the 
horizon limit would travel with the very outgoing null ray. Therefore, to obtain k[u) for 
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Figure 16. Temperature for the Unruh observer Kunruh a-s a function of r. We use 
2m = 1 units. 



the freely-falling observers from Kunmh, one has to extract from it a Doppler red-shift 
factor (going to zero at the horizon), or equivalently multiply it by a Doppler blue-shift 
factor (its inverse, diverging at the horizon). As we are going to show, this divergence 
counterbalance the zero value of ^unruh at the horizon yielding the appropriate finite 
result for the different khc's- 

Which is the precise form of this Doppler blue-shift factor? Our guess is 

d^ 



dr,.. 



( 

\dt 



dr,.„ 



(77) 



dt J \dt dt 

Here, r\{t) represents the trajectory of a null ray going away from the horizon and rrg(t) 
the freely-falling trajectory that starts with zero velocity at r = tq. In addition, we have 
put a subscript ro in ('^) to distinguish the different Doppler factors associated with 
freely-falling observers released at different radial distances rg. The reason behind this 
guess is the following. 

A physical Doppler factor has to compare the velocity of the massive objects with 
the velocity of light. 



D 



c — V 

C + V 



(7J 



However, note that using (t, r) coordinates the velocity of light is different from unity; 



that is why we have included dri/dt. The Doppler factor (77) can be alternatively 
written as 



1 - 



dr,,„ 



1/2 



1 + 



dr.^, 



-1/2 



(79) 



dri / \ dri 

In this way, one would compare the real local velocity of the observer with respect to 
the local velocity of light (which is always equal to one). 



We can take dr^o/dt from (56) and calculate dri/dt from the outgoing rays equation 
r*(ri) = const. 



dri 
d^ 
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(80) 
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After some manipulation, we finally find 
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If we compute now 

fiVnruh{r)Dro{r) 



(81) 



f82) 



and take the limit r — )■ 2m, we can check that one obtains precisely (71). So, as we 
advanced, it is really the extraction of a huge Doppler red-shift factor the reason why the 
freely-falling observers perceived some effective temperature even at horizon crossing. 
In other words, for long enough waiting times, a generic freely-falling observer sees a 
small but non-zero k when crossing the horizon. Only a hypothetical observer with zero 
radial velocity at the horizon (the Unruh observer there) would see a zero k, but because 
of the huge red-shift factor associated to its instantaneous trajectory (this observer is 
instantaneously travelling with the null ray). 

The final increase in n always present for freely-falling observers is caused by the 
Doppler shift factor. Take the case of an observer freely falling from infinity. As we 
explained, the perceived radiation starts at zero temperature. Then, if one waits long- 
enough, a plateau region is formed with an almost constant k. Finally, when the observer 
approaches the horizon, k increases up to the value 1/m, which corresponds to four times 
the standard Hawking temperature. The value of k from the appearance of the plateau 
on is due to the competition between Unruh's switching off of the radiation emission 
when going towards the horizon, which surprisingly passes through an intermediate 
increasing phase (see figure 16), and the multiplication by a Doppler blue-shift factor 
due to the velocity of the freely-falling observer with respect to the Unruh observer. 
The final result of this competition is the exact numerical computations we have shown 
throughout the paper (considering always sufficiently long waiting times). 

For instance, let us plot 

tiVnruh{r{u))Dro^ooir{u)) (83) 



against the exact numerical k,(u) of section 4.3 Here r{u) is the position of the observer 
falling from infinity as a function of the label u. From (22) one obtains r{u) 

Q n2/3 



r(u) = 2m 



Am 



(mh - n) + 1 



In figure 17, we compare this approximation with the exact result for k{u) for a long 
waiting time Ato = 200m. We can see that the approximation is perfect from the 
starting of the plateau till horizon crossing. Obviously, the previous approximation fails 
to capture the switching on of Hawking radiation. 

Actually, this perfect fit for long waiting times can also be checked analytically, 
and not only for observers falling from infinity, but for the generic case of freely-falling 



Hawking radiation as perceived by different observers 



28 




0.0 L 



-60 



-40 



-20 



Figure 17. Temperature as a function of u for a freely- falling observer from 

infinity with Atg — 200m using exact expression (curve depicted ) and using 

(Kunruh^ro^oo) ("") approximation (curve depicted ). We use 2ra — 1 units and 

mh = 0. 



observers from any radius. For instance, one just has to approximate the /€(?/, r) in (66) 
for Ato — c>o (this is equivalent to assuming U ~ Uh)-, which yields 
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and then see that it is precisely equal to 
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With hindsight, we could have used this equality to obtain all the previous partial 
results. But we can go now even further: if we put Kreieasedir, Ato) instead of Kunruh(^) 



in (86), then it is possible to see that we reproduce the entire formula for the effective 



temperature in (66). This is because between two freely-falling observers at the same 
position the only difference is a Doppler shift, and this fact has nothing to do with 
being in the Unruh state or not. We checked this first for the Unruh state because we 
wanted to explain the appearance of the final peak in k, a peak whose appearance was 
in apparent contradiction with the vanishing value we expected for that state. 
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6. Comparison with the Unruh state 



One might worry that some of the resuhs shown in this paper related with the late-time 
behaviour of the system, are peculiarities of an unusual choice of the vacuum state. 
However, this is not the case. In this section, we will briefly show that at late times 
this state is really equivalent to the Unruh vacuum state. Thus, all the properties found 
in the previous sections that applies for late times, are immediately present also in the 
Unruh state. 



f/n- So, as we did in (37), we can reverse 



The equivalence between the states for late times can be directly seen from (27) 
In that expression, late times means just U 
the approximate relation to find 

n 3 
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+ 1 



-ii/(4m) 
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But this is precisely the relation which defines the coordinate U that appears in the 
Unruh modes [H] (apart from the irrelevant origin of the coordinates). Thus, if we are 
using this coordinate to define the vacuum state, we are in fact choosing the Unruh 
state for late times. 

Nonetheless, just as an example, let us perform a specific calculation entirely based 
in the Unruh state choice. For instance, let us calculate the case of a freely-falling 
observer from infinity. With this aim, consider a redefinition of U using the usual and 
simple form 

U := -4me-^/(^'"\ (88) 



and take (44) as the definition of u. That is, U defines the vacuum and is the coordinate 



that appears in the Unruh modes; and u is the proper time that a freely-falling observer 



from infinity uses to label the rays. Replacing (44) in (88) we obtain a relation U{u) 
from which we can directly compute k,{u). The result is 



k(m) 
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From here one can check the limits 

1 



k{u — 7- 



-oo 



and /{(n — )■ mh) 



^9) 



(90) 



4m ~ ' m 

with the naked eye. The first limit is the usual Hawking radiation temperature, which 
in the Unruh state is always present, even at the asymptotic past. It is as if we extended 
the plateau to the u — )■ — oo limit (see figure 18). The second limit is equal to the one 



obtained at the horizon crossing for long enough delays in section 4.3 [see (52)], and 



also for far enough releasing radius in section 4.4 [see (72)]. Thus, the final peak in the 



temperature is not something coming from an eccentric selection of the state: it is a 
characteristic of the Unruh state. 
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Figure 18. Temperature as a function of u for a freely-falling observer from infinity 
in the Unruli state. We use 2m — 1 units and uh — 0. 

7. Summary and conclusions 

In this paper we have considered a quantum scalar field over a Schwarzschild black hole 
geometry. As a first step, we have set up the quantum scalar field in a specific vacuum 
state that interpolates between Boulware vacuum state at early times and Unruh vacuum 
state at late times. In other words, we have chosen a non-stationary vacuum state so 
that it simulates the physics of a dynamical collapse scenario. In this scenario, before 
the black hole forms there is no radiation at infinity but, after its formation, observers 
at infinity start receiving a Hawking fiux of radiation: the black hole Hawking emission 
switches on at some particular ignition time (of course, it is not an instantaneous event). 
This vacuum state has been fixed by requiring that a reference observer freely falling 
from infinity detects no radiation. 

As a second step, we have analysed the radiation associated with this vacuum 
state in terms of how it is perceived by different observers: static observers at a fixed 
radius; freely-falling observers from infinity, with different time delays with respect to 
the reference observer; and observers maintained at a fix radius for some time, and 
then, after being surpassed by the reference observer in its way in, released to fall freely 
towards the horizon at different releasing times. 

The method used follows the recent analysis in [U [13] . For each observer we define 
a function k varying along the observer's trajectory. We also define an adiabatic control 
function e as e = | d^/dti In the portions of the trajectory in which k stays almost 
constant (e <C 1) one can be sure that the observer detects thermal radiation at a slowly 
varying temperature /cbT = k/2'k. When this is not exactly the case, the radiation 
perceived is non-thermal, but still, we take the value of k as an estimate of the amount 
of perceived radiation and, somewhat abusing of the language, we talk about a time 
varying effective temperature. 

• Static observers at fixed radius. We show that they perceived an effective 
temperature that starts from zero, at early times, and after a transient period. 
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stabilizes at Hawking's temperature l/(87rm) multiplied by the gravitational blue- 
shift factor associated with the particular radial position of the observer. This 
description clearly reflects the ignition associated with the formation of a black hole 
we were seeking for. We also see that the transient period is quite sort: 0.07 x 10^'^ s 
for a collapsing neutron star with 1.5 times the mass of the Sun and a radius of 
12 km, when detected from far enough positions. Observers closer to the horizon 
will perceive this process with a slower pace due to the gravitational slow down of 
their clocks. 

• Observers freely falling from infinity. When their delay with respect to the reference 
trajectory used to fix the vacuum state is small, they perceive almost no radiation 
at all. The physical picture is the following: If an observer sitting at the very surface 
of a star that is collapsing in free fall to form a black hole perceived vacuum at the 
beginning of the collapse, he would perceive vacuum during his entire collapsing 
trajectory (this is reasonable, as he has not perceived any gravitational fleld). 

On the contrary, when the delay is long enough, the observer has time enough to 
see that the black hole is emitting Hawking radiation. He first sees the switching 
on of this radiation. Then, during some time he perceives an almost constant 
temperature radiation passing by. The temperature at this plateau is essentially 
the result of multiplying Hawking's temperature by a Doppler shift factor associated 
with the velocity of the observer. In this region the adiabatic condition is perfectly 
satisfied, so here to describe the perceived radiation as having a temperature is 
strictly correct. Finally, in the last stages of his approach to the horizon, and 
surprisingly at first sight, the effective temperature rises reaching exactly four times 
Hawking's temperature. The adiabatic condition is only mildly violated in these 
last stages, so although the perceived radiation is not exactly Planckian, it is not 
hugely different from Planckian either. 

• Observers freely falling from a radius. These observers are standing still at a fix 
radial distance ro from the horizon (supporting themselves against the gravitational 
pull with some rockets) till a release time (when the rockets are switched off), which 
happens some waiting time after the reference observer, in free fall, has overtaken 
them. If this waiting time is long enough, the observers will see Hawking radiation 
appear and settle down to a constant value, right in the same way as described 
above for static observers. Once the observer is released, the value of the effective 
temperature undergoes a jump with a value equal to the surface gravity associated 
with position ro times the gravitational blue-shift factor associated with this same 
position. We calculated the value of the effective temperature after the release and 
called it A;BTreieased('") = '«reieased('")/(27r). Prom the release time on, the effective 
temperature starts to rise, first slowly and later quicker, so as to reach a finite value 
at horizon crossing. The precise behaviour of the effective temperature depends on 
the release radius, and so does the asymptotic value of it, which becomes 4Th when 
To — >■ oo. 
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Once Hawking radiation is switched on, the value of feBTunmhl?") = /^Unruh('")/(2vr), 
that is, the effective temperature as perceived by an observer freely falling at r and with 
a vanishing instantaneous radial speed, interpolates from Th when r — )■ oo to zero when 
r — )■ 2m. This is consistent with the idea that the Unruh state is vacuum for freely 
falling observers at the horizon. However, in general one has to be careful with this 
assertion. We have shown that after Hawking radiation is switched on, generic freely- 
falling observers at the horizon do not perceive vacuum, but a finite effective temperature 
higher than their initial Unruh effective temperature. The reason is that the Unruh 
observer experiences a Doppler red-shift factor with respect to a generic observer (due 
to their velocity difference), or in other words, a generic observer experiences a Doppler 
blue-shift factor with the respect to the Unruh observer. This Doppler blue-shift factor 
diverges at the horizon in such a way that, when multiplied by the vanishing value of 
Unruh's temperature at the horizon, results in a finite temperature for a generic observer. 
In fact, we have proved the following completely general result: for any freely-falling 
observer the exact effective temperature T(r) perceived at position r is precisely the 
product of the temperature when released there Treieased('") by the Doppler blue-shift 
factor associated with the radial velocity of the observer when passing through that 
position. 

Regarding the observability of Hawking radiation, in pure theoretical terms the 
discussed results lead to the idea that the closer to the horizon one supports oneself 
against gravity, the better the chances of detecting its radiation, due to the gravitational 
blue-shift factor. However, in practical astrophysical scenarios the tiny Hawking 
temperature of a stellar mass black hole, about a few nano-Kelvins, will be masked by the 
~ 3 Kelvin of the cosmic microwave background. This radiation would also experience a 
gravitational blue shift if detected by observers close to the horizon, so that going closer 
to the horizon would not improve in principle the chances of detection. Another matter 
is the detection of Hawking radiation in analogue models of General Relativity. There 
exist quantum systems, as Bose-Einstein condensates, in which Hawking temperature is 
not that much masked by the environment temperature. In those situations one could 
think that analysing the presence of radiation close to the horizon would improve the 
detection chances, but this is not that clear. These systems exhibit modified dispersion 
relations so that when blue shifting the modes, they would enter into the non-relativistic 
regime, behaving very differently to what has been found and discussed in this paper. 
In addition, the size of the entire system in the lab will be so small that to control the 
physics localized close to horizon, a tiny fraction of the total size of the system, seems 
implausible. Other ways of detection that could prove better are the so called correlation 
measurements between outgoing and ingoing partners of the Hawking radiation [HI [15] . 

Let us finish the paper by mentioning that the results presented here will also 
constitute an adequate preparatory background material for a future work dealing with 
analyses of buoyancy over black hole horizons. 
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